----------------------------------------------------------------------------------
      name:  <unnamed>
       log:  G:\My Drive\Research\CTC expansion\submission\ntj\stata replication f
> iles/output/figure7.log
  log type:  text
 opened on:   3 Dec 2023, 18:09:18

. 
. version 16      

. 
. use "$data/acsmergeddata.dta", clear

. 
. sum newrtwchange if calckids>0 [aw=perwt]

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange | 356,390    39120804    .0408424   .0405399          0   .2283467

. bysort married: sum newrtwchange if calckids>0 [aw=perwt]

----------------------------------------------------------------------------------
-> married = 0

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange | 122,494    15347076    .0637255   .0502828          0   .2283467

----------------------------------------------------------------------------------
-> married = 1

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange | 233,896    23773728    .0260703    .022716          0     .21453


. 
. sum newrtwchange if calckids>0 & sex==2 & agi<80000 & incwage>0 [aw=perwt]

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange | 182,692    20319255    .0546366   .0451657          0   .2283467

. bysort married: sum newrtwchange if calckids>0 & sex==2 & agi<80000 & incwage>0 
> [aw=perwt]

----------------------------------------------------------------------------------
-> married = 0

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange |  70,949     8956824    .0870148   .0432268          0   .2283467

----------------------------------------------------------------------------------
-> married = 1

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange | 111,743    11362431    .0291135    .026405          0     .21453


. 
. sum newrtwchange if calckids>0 & sex==2 & agi<40000 & incwage>0 [aw=perwt]

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange | 145,354    16185245    .0563398   .0476621          0   .2038863

. bysort married: sum newrtwchange if calckids>0 & sex==2 & agi<40000 & incwage>0 
> [aw=perwt]

----------------------------------------------------------------------------------
-> married = 0

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange |  53,053     6804768     .095872   .0415176          0   .2038863

----------------------------------------------------------------------------------
-> married = 1

    Variable |     Obs      Weight        Mean   Std. dev.       Min        Max
-------------+-----------------------------------------------------------------
newrtwchange |  92,301     9380477    .0276623    .026704          0     .19706


. 
. replace newrtwchange=max(-20,-100*newrtwchange)
(323,834 real changes made)

. 
. twoway hist newrtwchange [fw=perwt] if calckids>0 & married, color(gs6) w(1) fra
> c ||  hist newrtwchange [fw=perwt] if calckids>0 & married==0, w(1) frac fcolor(
> none) lcolor(black) graphregion(color(white)) bgcolor(white) ytitle("Fraction of
>  Sample") text(.3 -10 "Average = -4.1%" "Married Average = -2.6%" "Unmarried Ave
> rage = -6.4%" , j(left) size(medsmall)) xtitle("Change in Return to Work (Pct Pt
> s)" "{bf:Panel A: All Parents}") xlabel(0(-4)-20) ylabel(0(.1).35) leg(order(1 "
> Married" 2 "Unmarried") c(1) pos(11) ring(0))

. graph save "$output/newrtwchangechangehist1.gph", replace
(file G:\My Drive\Research\CTC expansion\submission\ntj\stata replication
    files/output/newrtwchangechangehist1.gph not found)
file G:\My Drive\Research\CTC expansion\submission\ntj\stata replication files/out
> put/newrtwchangechangehist1.gph saved

. 
. twoway hist newrtwchange [fw=perwt] if calckids>0 & married & agi<80000 & incwag
> e>0, color(gs6) w(1) frac ||  hist newrtwchange [fw=perwt] if calckids>0 & marri
> ed==0 & agi<80000 & incwage>0, w(1) frac fcolor(none) lcolor(black) graphregion(
> color(white)) bgcolor(white) ytitle("Fraction of Sample") text(.3 -10 "Average =
>  -5.5%" "Married Average = -2.9%" "Unmarried Average = -8.7%" , j(left) size(med
> small)) xtitle("Change in Return to Work (Pct Pts)" "{bf:Panel B: Working with A
> GI < $80,000}") xlabel(0(-4)-20) ylabel(0(.1).35) leg(order(1 "Married" 2 "Unmar
> ried") c(1) pos(11) ring(0))

. graph save "$output/newrtwchangechangehist2.gph", replace
(file G:\My Drive\Research\CTC expansion\submission\ntj\stata replication
    files/output/newrtwchangechangehist2.gph not found)
file G:\My Drive\Research\CTC expansion\submission\ntj\stata replication files/out
> put/newrtwchangechangehist2.gph saved

. 
. twoway hist newrtwchange [fw=perwt] if calckids>0 & married & agi<40000 & incwag
> e>0, color(gs6) w(1) frac ||  hist newrtwchange [fw=perwt] if calckids>0 & marri
> ed==0 & agi<40000 & incwage>0, w(1) frac fcolor(none) lcolor(black) graphregion(
> color(white)) bgcolor(white) ytitle("Fraction of Sample") text(.3 -10 "Average =
>  -5.6%" "Married Average = -2.8%" "Unmarried Average = -9.6%" , j(left) size(med
> small)) xtitle("Change in Return to Work (Pct Pts)" "{bf:Panel C: Working with A
> GI < $40,000}") xlabel(0(-4)-20) ylabel(0(.1).35) leg(order(1 "Married" 2 "Unmar
> ried") c(1) pos(11) ring(0))

. graph save "$output/newrtwchangechangehist3.gph", replace
(file G:\My Drive\Research\CTC expansion\submission\ntj\stata replication
    files/output/newrtwchangechangehist3.gph not found)
file G:\My Drive\Research\CTC expansion\submission\ntj\stata replication files/out
> put/newrtwchangechangehist3.gph saved

. 
. graph combine "$output/newrtwchangechangehist1.gph" "$output/newrtwchangechangeh
> ist2.gph" "$output/newrtwchangechangehist3.gph", col(1) graphregion(color(white)
> )

. graph export "$output/figure7.png", replace
file G:\My Drive\Research\CTC expansion\submission\ntj\stata replication
    files/output/figure7.png saved as PNG format

. 
. erase "$output/newrtwchangechangehist1.gph"

. erase "$output/newrtwchangechangehist2.gph"

. erase "$output/newrtwchangechangehist3.gph"

. 
. log close
      name:  <unnamed>
       log:  G:\My Drive\Research\CTC expansion\submission\ntj\stata replication f
> iles/output/figure7.log
  log type:  text
 closed on:   3 Dec 2023, 18:09:24
----------------------------------------------------------------------------------
